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ABSTRACT 

Context. The ionizing Lyman continuum flux escaping from high redshift galaxies into the intergalactic medium is a fundamental 
quantity to understand the physical processes involved in the reionization epoch. However, from an observational point of view, 
direct detections of HI ionizing photons at high redshifts are feasible for galaxies mainly in the interval 2^3 — 4. 

Aims. We have investigated a sample of star-forming galaxies at 2 ~ 3.3 in order to search for possible detections of Lyman 
continuum ionizing photons escaping from galaxy halos. 

Methods. UV deep imaging in the COSMOS held obtained with the prime focus camera LBC at the LET telescope was used 
together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 
45 galaxies at 2 ~ 3.3 with L > 0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 
3.27 < 2 < 3.40 both in the R and deep U bands (magnitude limit U 29.7 at S/N=l). At these redshifts the R band samples the 
non-ionizing 1500 A rest frame luminosity and the U band samples the rest-frame spectral region just short-ward of the Lyman edge 
at 912 A. Their flux ratio is related to the ionizing escape fraction after statistical removal of the absorption by the intergalactic 
medium along the line of sight. 

Results. A sub-sample of 10 galaxies apparently shows escape fractions > 28% but a detailed analysis of their properties reveals that, 
with the exception of two marginal detections {S/N 2) in the U band, all the other eight galaxies are most likely contaminated 
by the UV flux of low redshift interlopers located close (in angular position) to the high-z targets. The average escape fraction 
derived from the stacking of the cleaned sample was constrained to fff) < 2%. The implied hydrogen photo-ionization rate is a 
factor two lower than that needed to keep the intergalactic medium ionized at 2 ~ 3, as observed in the Lyman-a forest of high-z 
QSO spectra or by the proximity effect. 

Conclusions. These results support a scenario where high redshift, relatively bright {L > 0.51/*) star-forming galaxies alone are 
unable to sustain the level of ionization observed in the cosmic intergalactic medium at 2 ~ 3. Star-forming galaxies at higher 
redshift and at fainter luminosities {L « L*) can be the major contributors to the reionization of the Universe only if their 
physical properties are subject to rapid changes from 2 3 to 2 ~ 6 — 10. Alternatively, ionizing sources could be discovered 

looking for fainter sources among the AGN population at high redshift. 

Key words. Galaxies: distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: photometry 


1. Introduction 

The reionization epoch marks a fundamental event in 
the thermal history of the Universe, probably located 
in the redshift interval 2 = 6 — 10. The lower limit 
is derived from observations of the Gunn-Peterson effect 
in z=6 QSO spectra (IFan et al. 2006|l while the upper 
end comes from measurements of the Thomson scatter- 


Send offprint requests to: A. Grazian, e-mail: 

andrea.grazian@oa-roma.inaf.it 
* Based on data obtained with the European Southern 
Observatory Very Large Telescope, Paranal, Chile, under Large 
Program 185.A-0791 and on observations made at the Large 
Binocular Telescope (LBT) at Mt. Graham (Arizona, USA). 


ing in the CMB polarization map by Wilkinson Microwave 
Anisotropy Probe (WMAP, IHinshaw et al. 201311 and 
Planck (|Ade et al. 2ni5|l . The time extension of the reion¬ 
ization processes can also be constrained by the so-called 
kinetic Sunyaev-Zel’dovich (kSZ) effect. Combining data 
from the South Pole Telescope (SPT) with the WMAP re¬ 
sults, Zahn et al. (2012) provided a 95% c.l. limit to the 
duration of reionization of Az < 4.4, placing the end of 
reionization at z = 7.2. While we are able to approxi¬ 
mately locate the reionization epoch, we still do not have a 
clear perception of the sources providing the bulk of the UV 
photons required to reionize the Universe. Prime suspects 
have been searched so far among high redshift star forming 
galaxies or AGNs. 
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It is well known that bright QSOs are effi¬ 
cient producers of HI ionizing photons (A < 912 

A), and they can ionize large bubbles of hydro¬ 
gen and helium even at dis tances up to several 
Mpcs (Cowie et al. (2009)[ Prochaska et al. (2009)[ 


same area of Shapley et al. (2006) [ namely the SSA22 field. 


Worseck et al. (2014a)| |Stevans et al. (2014)p. However, 

their volume density is apparently too low at high redshifts 
(z > 2) to provide the cosmic photo-ionization rate re¬ 
quired to keep the intergalactic medium (IGM) ionized at 
z ^ 6 (|Fan et al. 2006|l . Recent searches for low luminosity 
AGNs at high redshifts by means of deep X-ray imaging 


They found a significant population of LyC emitters at 
z > 3, with escape fractions in the range ~ 4 — 100% 
(some Lyman-a emitters with > 100%), providing the 
first indication that the escape fraction of galaxies started 
to increase at z > 3 (or that their sample may be con¬ 
taminated by foregr ound objects, as we will discuss in the 


following). Recently, Mostardi et al. (2013) studied a new 
sky region around an overdensity of sources at z ~ 2.8, 
reaching slightly similar conclusions {fesc ~ 5 — 49%). 
Results from other teams based on very deep nar- 


Giallongo et al. (2012) 

Giallongo et al. (2015)P leaving per limits of /J®^ < 5 - 10% to 

the escape fraction of 

Vanzella et al. (2010b) 

room for important contributions to the reionization bv star tormina galaxies at 2 = 3 — 4 I 


this fainter AGN population. 

High redshift star forming galaxies can also be im- 
portant contributors of ionizing photons at z > 2 

( Haardt fc Madau (2012)[ [Duncan fc Conselice 2015|l . An 
intense star formation activity is observed till z 8 — 10 


see for example McLeod et al. (2015)1, and the luminos¬ 


ity function (LF) of galaxies at z > 4 is gra dually steep- 


percents or less (Grimes et al. (2009)| 

Cowie et al. (2009) 

Siana et al. (2010)1 

Bridge et al. (2010) 

Leitet et al. (2013)1 

. At high redshifts the current sit- 


star-forming galaxies measured by different groups in the 
last 15 years led to different results indicating that some 
systematics are still present in the analysis and in the 
galaxy samples used. 


ters by Iwata et al. (2009) and Nestor et al. (2011) The 


original result of|[Shapley et al. (2006) should thus be con- 


verted into a strict upper limit of < 2% at z ~ 3 


(Cooke et al. (2014)I. 

With the aid of n arr ow band imaging at A ^ 3590 
A, Iwata et al. (2009) and Nestor et al. (2011) studied the 


Boutsia et al. (2011)). The main reason for the discrepant 


ening with redshift (Finkelstein et al. (2014) I. Young stel¬ 
lar populations inside high redshift galaxies are produc¬ 
ing a large amount of UV photons, but only a fraction of 
them can reach the surrounding IGM. An efficient way of 
parametrizing the emission of HI ionizing photons w.r.t. 
the non-ionizing ones is the definition of the escape fraction 


(Steidel et al. (2001)). The escape fraction parameter is a 
measurement of the fraction of HI ionizing photons which 
are not absorbed by the interstellar medium (ISM) and are 
thus free to ionize the neutral hydrogen in the IGM. The 
extensive definition of the escape fraction will be provided 
in Sect. 3. 

The bright star forming galaxies at low redshifts 
(z < 2) in general are not able to ionize their neigh¬ 
borhoods since their escape fractions are only few 


results is due to the details of the photometric analysis 
and to the different morphological analysis of the candi¬ 
date Lyman continuum emitters. 

The most convincin g explanation for this d iscrepancy 
has been proposed by Vanzella et al. (2010a) (hereafter 
VlOa) on the basis of possible UV hux contamination in 
the Nestor’s and Mostardi’s LyG emitting sample by low 
redshift interlopers which are near the position of the high 
redshift galaxies. As estimated by VlOa, the chance for con¬ 
tamination increases with large seeing and deep UV images. 
Gonsidering a sample of 101 galaxies at 3.4 < z < 4.4 in 
the GOODS-South field they found only one convincing de¬ 
tection after removing probable contaminants identified on 
the basis of accurate multicolor HST photometry at high 
angular resolution. 


At the same time Boutsia et al. (2011) studied a small 
sample of 11 galaxies at z ~ 3.3 in the Q0933, Q1623, and 
GOSMOS fields with very deep imaging in the U and R 
bands obtained at the prime focus camera LBC of LBT. 
The non detection of all the galaxies resulted in a strict 
upper limit to the escape fraction of < 5%, but the 
result was affected by low number statistics and by the un¬ 
certainty due to the stochasticity of the IGM absorption. 
Their conclusions were that, at present, a class of ioniz¬ 
ing sources responsible for cosmic reionization has not been 
clearly identified and a step forward to the escape fraction 
studies is required. _ 


Steidel et al. (2001) claimed the detection of significant 
emission (Z' ®^ ^ 507o) of Lyman Gontinuum (LyG) photons 
from the composite spectrum of 29 Lyman break galax¬ 
ies (LBG) at z ^ 3. Two LBGs from the Steidel high- 
redshift sample wer e spectroscopically obser ved at interme¬ 
diate resolution by Giallongo et al. (2002)| with very deep 
integration. They derived only an upper limit to the rela¬ 
tive escape fraction of < 16%. Further spectroscopic 
work by Steidel’s team found only two emitting galaxies at 
z ~ 3, out of 14 analyzed, with significant escape fraction 
{fisc > 50 — 100%). For the remaining population only up¬ 
per limits were derive d resulting in an avera ge /J®^ ~ 14% 
for the whole sample (Shapley et al. (2006)1. However, the 
two emitting galaxies were not conhrmed as LyG emit- 


Recently, Siana et al. (2015) re-analyzed with the HST 
instrument WFG3-UVIS 5 galaxies that were identified 

LyC candidates at z ~ 3.1. 
images and NIR Keck spec- 


as 


by Nestor et al. (2011) 

Thanks to deep 
tra of these galaxies, Siana et al. (2015) concluded that 
the supposed LyC emission is in fact provided by lower 
redshift interlopers. The new upper limit by these ob- 
servati ons is /J®^ 7 — 9% , con sistent with the lim- 

its by Vanzella et al. ( 2010b)| and Boutsia et al. (2011) 


In addition, [Mostardi et al. (^015) reported an HS'L fol 


low up analysis of 16 candidate LyC emitters found by 


Mostardi et al. (2013) 
Jl25 


They obtained deep C/ 336 , Veoej 
and H 160 imaging with WFC3 and ACS for these 


galaxies: thanks to the high spatial resolution of HST, they 
showed that 15 emitters can be explained with contami¬ 
nation of low-z objects, wrong spectroscopic redshifts or 
uncertain photometric redshifts. Only one galaxy shows a 


robust LyC detection, according to Mostardi et al. (2015) 
On the basis of these results, they concluded that galax- 
ies at z ^ 3 are providing the same contribution of QSOs 
to the HI photo-ionizing background, although with large 
uncertainties. 
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Increasing the number of accurate LyC measures is not 
a trivial task since it requires either very deep spectroscopy 
or very deep UV imaging of a large sample of galaxies 
with known spectroscopic redshifts in the range z = 3 — 4 
where the IGM still allows t he transmission of so me ionizing 
flux along the line of sight (|Inoue et al. (2014)|). Indeed at 
higher redshifts it is difficult to disentangle the absorption 
due to the galaxy ISM from that due to the surrounding 
IGM along the line of sight (LoS). 

In this work we exploit the availability of very deep U- 
band photometry obtained with the LEG prime focus im¬ 
ager at LET, coupled with a large number of galaxies with 
spectroscopic redshifts in a relatively narrow range around 
z ~ 3.3 that allows the UV band to sample the spectral re¬ 
gion just short-ward of 912 A rest. The spectroscopic sam- 
ple is taken from the VUDS survey (Le Fevre et al. (2015)) 
of high redshift galaxies (2 < z <~ 6) in the COSMOS field 
where deep HST imaging is also available and allows to look 
for possible contamination by low redshift interlopers. 

This paper is organized as follows. In Sect. 2 we present 
the dataset, in Sect. 3 we describe the method adopted, in 
Sect. 4 and 5 we show the results for individual objects and 
for the overall sample as a whole, in Sect. 6 we discuss our 
results and in Sect. 7 we provide a summary. Throughout 
the paper we adopt the A-CDM concordance cosmological 
model (Ho = 70 km/s/Mpc, VLm = 0.3 and Ua = 0.7). All 
magnitudes are in the AE system. 


2. Data 


(Grogin et al. (2011)1 
imaging dataset of tin 


in detail in Eoutsia et al. (2011) and Eoutsia et al. (2014) 
The LEG instrument (Giallongo et al. (2008)1 has two 


(Eoutsia et al. (2011)1, using the same technique that will 
be adopted in this work. 

The U-band observations of the GOSMOS field by LEG 
are constituted by three partially overlapping pointings of 
^ 6 hours each, for a total of 11 hours of exposure on the 
deepest regions. The relative good seeing (^ 1.0 arcsec) and 
the high sensitivity of the LEC camera in the UV allow us 
to reach a depth of 28.7 AE magnitudes (at S/N=l) over an 
area corresponding to 700 sq. arcmin., with slightly deeper 
limits ('^ 29.1) on the regions covered twice ('^ 450 sq. 
arcmin). 

The R band dataset has an average seeing of 0.9 arc- 
sec and in 3.3 hours of exposure time it reaches, for two 
third of the area covered by deep U band imaging, a depth 
of R=25.7 AE magnitude at S/N=10. The third pointing 
is shallower (R=23.8), since the exposure time was shorter 
(450 seconds). In this sub-region, we select only galaxies 
down to i? < 24.0 to work at a good S/N ratio in the detec¬ 
tion band, which will be used to measure the non ionizing 
continuum of the sources. 

has been used to 

__ _ described in de- 

Giallongo et al. (2008) [ Eoutsia et al. (2011) and 
In particular, after correcting the 
images lor bias subtraction, we have applied standard 
flat-fielding using sky flats obtained during twilight. 
The sky background was estimated using SExtractor 
([Eertin fc Arnouts 1996|) with a mesh of 64 by 64 pixels 
and adopting a median of 3 by 3 meshes. The astrometric 
corrections were computed using the software “AstromC” 


The standard LEC pipeline 
reduce the imaging dataset, as 


tail in _ 

Eoutsia'et al. (2014) 


The COSMOS field (Scoville et al. (2007)1 has been se¬ 
lected as an optimal area for the measurement of the es¬ 
cape fraction of galaxies at z ^ 3, due to its combination 
of deep multiwavelength imaging by LEC coupled with a 
large number of spectroscopic surveys over a broad redshift 
interval 0 < z < 7 obtained by extensive observat ional 
campaigns ( Lilly et al (2007)] Le Fevre et al. (2015) I. 

In the present anal ysis we have also includ ed the galax¬ 
ies already studied by Eoutsia et al. (2011) in the Q0933 
and Q1623 helds. 

2.1. The photometric dataset 

The imaging data set used in this paper consists of deep 
LEC observations of the COSMOS field in the U and R fil¬ 
ters. These deep observations are part of a more ambitious 
program, the LEC/CANDELS ( [Grazian et al. in prep.[ ) 
survey, with the goal of providing complementary 
deep UV coverage with LEC to the northern fields 
(COSMOS, GOODS-North, EGS) of the CAND ELS survey 


described by Radovich et al. (2004) After resampling all 


the individual images in a common reference frame, we co¬ 
added the various images to create deep mosaics for each 
field. 

For each LEC image, we have derived a careful estimate 
of the absolute rms map taking into account the Poisson 
statistics and the instrumental gain for each individual chip. 


as described in Eoutsia et al. (2011) We have then prop¬ 


agated the rms map over the whole data reduction chain, 
obtaining as a final product of the pipeline an absolute vari¬ 
ance map associated with each stacked scientific frame. 


2.2. Spectroscopic redshifts 


Koekemoer et al. (2011) I. The LEC 
COSMOS held has been described 


wide-held imagers, one optimized for the U V wavelengths 
and the other for the optical-red bands. This imager 
has observed part (^ 0.20 sq. deg.) of the COSMOS 


(Scoville et al. (2007)1 held in UGRIZ bands to search for 
LEGs at 2.7 < z < 3.4 ([Eoutsia et al. (2014)|). These data, 
obtained with one of the most efficient ground-based large 
held UV imager at an 8 meter class telescope, the LET, 
have allowed us to derive a stringent limit of < 5% to 
the escape fraction of 11 star-forming galaxies at z ^ 3.3 


The VUDS spectroscopic survey (Le Fevre et al. (2015) I is 
an ESO Large Program (P.I. Le Fevre) aiming at obtain- 
ing ^10,000 spectroscopic redshifts with the VIMOS instru¬ 
ment to study early phases of galaxy evolution in the red¬ 
shift range 2 < z < 6 and beyond. This ambitious program 
has been conducted over 3 helds (VVDS-02h, COSMOS, 
ECDFS) and the targets were selected with different cri¬ 
teria (inclusive combination of U-dropout, E-dropout, and 
photometric redshift selections), tailored to maximize the 
efficiency of hnding galaxies at 3 < z < 5, thus being com - 
plementary to the zCOSMOS survey (Lilly et al. (2007)). 
The spectra are the result of deep integration (14 hours per 
target) covering the wavelength range 3650 < A < 9350 A. 
Details on the observational strategy, candidate selection, 
data reduction, quality hag assignment and hrst results of 


the VUDS survey can be found in Le Fevre et al. (2015) 
The data acquisition, reduction and redshift extraction pro- 
cedure of the whole VUDS survey have been completed re¬ 
cently. 
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We focus on a very small redshift range (3.27 < 
< 3.40) in order to minimize the effect of the strong 


IGM opacity evolution with redshift (Inoue et al. (2014) 
Worseck et al. (2014a) I and allowing a direct measure¬ 
ment of the relative ionizing escape fraction from star 
forming galaxies at 2 = 3.3. As explained also in 

Boutsia et al. (2011)| this narrow redshift range allows us 
to use the LBC U-band filter (SDT-Uspec) to sample the 
rest frame wavelengths short-ward of 912 A, but not too 
far from it, namely at > 880 A rest frame. An extension of 
the redshift range at z > 3.40 is possible, but it requires 
detail simulations of the variance due to the stochasticity 
of intervening dampe d Lyman-g systems or Lyman Limit 

Since it goes 
le analysis to 


filter of LBC and ado pting the same procedure outlined in 
Boutsia et al. (2011) we obtain an IGM mean transmission 
< exp{—Tg^^^} >= 0.28 at z ~ 3.3, similar to that foun d 


adopting the IGM corrections by Prochaska et al. (2009) 

It is worth noting that we are not considering here the 
variance of the IGM absorption in different lines of sight due 
to a patchy pattern of LLS. Our relatively large number of 
galaxies allows us to deal, from the statistical point of view, 
with the possibility that few lines of sight are completely 
absorbed by intervening optically thick systems at z < 3.3. 
In the discuss ion section we will q uanti fy this probability. 
Following Steidel et al. (200T)] and Siana et al. (2007) 


systems, as shown in Vanzella et al. (2010b) 
beyond the scope of this paper, we limit t 
z < 3.40. 

We started from the sample of 73 galaxies with highly 
reliable (flag=3, 4, 23, 24) spectroscopic redshifts within 
the interval 3.27 < z < 3.40 in the COSMOS area cov¬ 
ered by the VUDS spectroscopic survey. We do not con¬ 
sider within our sample the AGNs from the optical clas¬ 
sification of the VUDS spectra (flag=ll, 12, 13, 14, 19 or 
211, 212, 213, 214, 219 due mainly to the presence of a 
strong CIV emission line). We excluded also the spectro¬ 
scopic redshifts based only on a single emission line (flag=9 
or 29). Of the pure (non AGN) galaxies, only 35 are covered 
by the LBC images in the U and R bands, which do not 
cover the whole COSMOS field, as mentioned above. One 
of these sources has been excluded since it falls on an hori¬ 
zontal stripe of bad pixels caused by a saturated star in the 
LBC ima ges. Eleven additional galaxies were already dis¬ 
cussed by Boutsia et al. (2011)| resulting in a total sample 
of 45 star forming galaxies at z ^ 3.3 with secure spectro¬ 
scopic redshifts and deep U and R band images by LBC. 
Table [1] summarizes the properties of the individual galax¬ 
ies which are used in this paper to study the escape fraction 
at z = 3.3. 


3. The adopted method to derive the LyC escape 
fraction of star-forming galaxies in the 
VUDS-LBC/COSMOS field 

In this paper we exploit the deep UV imaging by LBC in the 
VUDS/COSMOS field to derive a meaningful measurement 
of the ionizing radiation of star forming galaxies at z ~ 3.3. 

The LyC escape fraction of high-z galaxies can be eval¬ 
uated by adopting very deep UV imaging data sampling 
the A < 912 A rest frame (Lyman continuum) region. Since 
the intrinsic spectral energy distribution (SED) of a galaxy 
is not known a priori, especially in the rest-frame far UV 
where significant dust reddening could be present, a rela- 
tive escape fraction is usually derived from the observations 
(Steidel et al. (2001)). It is defined as the fraction of escap- 
ing Lyman continuum photons divide d by the fraction o f 
escaping photons at 1500 A rest frame (Siana et al. (2007) I. 

The SED of a high-redshift source at A < 912 A rest 
frame is also affected by the photoelectric absorpt ion by 
the IGM (Moller & Jakobsen (1990) Madau (1995) I and it 
is mainly due to Lyman Limit Systems (LLS) with column 
densities of neutral hydrogen between logi^Nni) ~ 17 and 
20. To correct for the mean effect of the IGM absorption 


we estimate the relative escape fraction using the following 
equation: 


j;rel {Ll^mIL'im)int t_IGM\ 

/esc— !-TJ /TP \ ^‘'^Pvgoo > 


{Fii/Fu)oba 

which is equivalent to: 


( 1 ) 


= exp(-Tjj7)xl0-°-4(^™‘>-^i=“) = //j'/xl0°'^^i"™(2) 

where (Lisoo/Agoo)™* is the ratio of the intrinsic lumi¬ 
nosities at 1500 and 900 A rest frame, {Fn/Fu)obs is the 
ratio of the observed fluxes in the R and U band, exp(—r^f/) 
indicates the ISM extinction, exp(—r^f/) x is 

the absolute escape fraction /“s/, A 1500 and Agoo are the 
dust absorpt ion coefficients at 1500 and 900 A, respectively. 
According to|Calzetti et al. (2000)1 ^1500 = 10.33if(i? —U) 
and Agoo = 18.02A'(i^ - l^), thus giving = 2.59 for 
E{B — V) =0.1. According to Eq. 2, this implies that 
is always less than 100%. 

We adopt an intrinsic ratio of (Li 5 oo/A 9 oo)mt = 3 for 
our star-forming galaxies, taking into account the absorp- 
tion produced by th e st ellar atmospheres, as adopted by 
Steidel et al. (2001)| and |Boutsia et al. (2011)| on the ba- 
sis of predictions by the spectral synthesis stellar models. 
The intrinsic ratio (Lisoo/Agoo)™* depends on the physical 
properties of the galaxies, mainly the mean stellar ages, 
metallicities, IMFs, and SFHs. For typical star forming 
galaxies at z ~ 3.3 the intrinsic ratio (Li 5 oo/A 9 oo)mt varies 
between 1.7 and 7.1 for age between 1 Myr and 0.2 Gyr, 
adopting the BC03 library (IBruzual fc Chariot 2003| l. In 
the following sections we will discuss the implications of 
our choice, and we will show that our final results on the 
ionization emissivity are robust against the adopted value 
of 3 for the intrinsic ratio (Li 5 oo/Agoo)mt- To derive the 
observed fluxes at 900 and 1500 A rest frame at z ~ 3.3 
we rely on the U and R band flux ratio measured by LBC, 
{Fji/Fu)obs- 


we adopt the recent estimates by Worseck et al. (2014a) 
Convolving the IGM absorption at z = 3.3 with the U-banc 


In our previous paper (Boutsia et al. (2011)), we have 
estimated fluxes using aperture photometry matched to the 
seeing of the U and R bands. Here we explore a slightly 
different approach: if LyC photons are emitted by a par¬ 
ticular region of a galaxy characterized by low ISM ab¬ 
sorption, then we must observe also non-ionizing emission 
emerging from the sam e spatial region (see section 4.1.1 of 
Vanzella et al. (2012a) hereafter V12). As a consequence, 
the emission in the U band (sampling ~ 900 A at z ~ 3.3) 
must be always associated with an emission in the R band 
(~ 1500 A at z ~ 3.3). The opposite would not be true, 
however, since non-ionizing radiation can be detected in 
regions of a galaxy where the ISM or dust absorption are 
elevated, with a negligible emission of LyC photons. The 
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observed profile of a galaxy in the R band can thus be 
used to measure the ionizing radiations in the U band 
only in the physical regions where such photons are ex¬ 
pected, limiting the contamination by foreground galaxies 
lying by chance on a given line of sight of a star forming 
galaxy at z ^ 3.3. We t hus adopt the ConvPhot software 
(De Santis et al. (2007)1 to derive less contaminated pho- 
tometry in the U band, using the R-band image, with its 
much higher S/N ratio, as a prior to reconstruct the profile 
of galaxies in the U band. Although the computed magni¬ 
tude is not significantly different from that derived by aper¬ 
ture photometry, the photometric uncertainties are smaller 
in the ConvPhot procedures since the flux errors in the U 
band are, at least at a first approximation, equal to the re¬ 
sult of the sum over its rms^ weighing for the profile of the 
obje ct in the R band. Com paring the objects in common 
with Boutsia et al. (201 1)[ we find that this new method 
allows a gain by a factor of ^1.7 in S/N ratio in the U 
band photometry, and thus a similar improvement on the 
determination of the uncertainties on the LyC escape frac¬ 
tion. The resulting limiting magnitude in the LBC U band 
image is thus ~ 29.7 AB at S/N=l for the central overlap¬ 
ping region, and slightly shallower (^ 29.3) in the flanking 
areas covered with only one LBC pointing. 


Following V12, we define “local” escape fraction as the 
parameter measuring the LyC emission w.r.t. non ionizing 
radiation, measured in the same physical region where both 
the 900 and 1500 A rest frame wavelengths are emitted. As 
a sanity check for the consistency of the following results, 
the “local” escape fraction cannot exceed 100% for reason¬ 
able values of the IGM transmission and of the intrinsic 
ratio (Li 5 oo/L 9 oo)mt- As in V12, a “local” escape fraction 
of 500-1000% can be a clear indication of contamination in 
the LyC band from a foreground interloper. 


4. The relative escape fraction of individual 
galaxies at Z/ 3.3 


Using the method outlined in the previous section, we com¬ 
puted the relative escape fraction of the 45 individual galax¬ 
ies at 3.27 < Zspec < 3.40 of the VUDS-LBC/COSMOS 
sample. In case of non detection in the U band (i.e. SjN < 
1 ) we provide the upper limit of the relative escape fraction 
at Icr. This is a fair limit since we are forcing the measure¬ 
ment of LvC flux in a region where we see radiation at 
A ^ 1500 A rest frame (non ionizing) from the clear detec¬ 
tion in the R band. Table [I] summarizes the properties of 
these objects. As we will discuss in section 6.2, the patchy 
nature of the IGM makes the determination of the escape 
fraction for individual galaxies slightly uncertain. Our main 
results has thus been derived in section 5.1 by studying the 
properties of the whole galaxy population at 2 ^ 3 through 
image stacking. 


Table[T]shows that there are 10 detections in the U band 
within the initial sample of 45 star forming galaxies at z ~ 
3.3 in the VUDS area covered by L BC observations (com - 
plemented with early results from Boutsia et al. (2011) I. 
At face value, this result suggests possible LyC emission 
for this sub-sample. 


4.1. Galaxies with possible LyC detection 

We discuss here the properties of the ten interesting objects 
at z ~ 3.3 which are detected in the U band of LBC showing 
possible ionizing emission. In the next sub-sections we check 
the reliability of the LyC sources based on the deep HST 
F814W (hereafter indicated as I814W) band and the LBC 
U and R bands. 


4.1.1. Galaxies with //j/ > 100% 


In Table[T]there are four galaxies with a non-physical escape 
fraction ranging between 140 and 292% (last four objects 
at the end of the table). 

Fig[T] shows as an example the cutouts in the LBC U 
and R bands, and the HST image in the I814W filter of 
the galaxy with ID=16669 in Table [TJ The spectroscopic 
redshift from VUDS is associated with a compact galaxy 
(indicated by the green circle in the HST image on the 
right). The compact source falls near (distance=2.1 arcsec) 
a bright and extended object which is clearly detected in the 
U band, while at the position of the z 3.3 galaxy no flux 
is emitted in the U band. When computing the escape frac¬ 
tion with our method (i.e. detecting the galaxy in the LBC 
R-band and measuring the LyC flux in the LBC U-band on 
the same area, shown by the red ellipse in the central panel 
of FiglD), the two galaxies are merged in a single detection 
in the R band, and thus our procedure gives a significant 
LyC detection from the U band. Translating the observed 
U and R fluxes into an escape fraction, it turns out that 
this object has = 229.5%, more than twice the maxi¬ 
mum value allowed for a typical star forming galaxy. If we 
restrict the investigated area to the region with a significant 
U-band emission (the area of the red ellipse in FiglHoutside 
the green circle; this barely extended object is also visible 
in the HST image), and we compute the flux ratio between 
the U and R band, we end up with an escape fraction of 
519.4%. The uncertainties in the U and R magnitudes are 
of the order of ^ 0.05, indicating that the latter estimate of 
the “local” escape fraction is robust. The extended galaxy 
on top of the compact z 3.3 source is 3 arcsec in diame¬ 
ter, which corresponds t o 22 kpc if the two gala xies are at 
the same redshift, while Bouwens et al. (2004) found that 
a typical L* star forming galaxy has an half light radius of 
2 kpc at z ^ 2.5. A plausible explanation for the detection 
in LyC for source ID=16669 is thus the presence of a fore¬ 
ground object which mimics the presence of a LyC region 
near the z ^ 3.3 galaxy. 

A similar situation happens for sources ID=52336 and 
55594 lFig lA.il) . where the presence of an extended galaxy 
and a compact object, respectively, with relatively bright 
U band emission but with significant offset (distance of 
2.4 and 1.5 arcsec, respectively) from the z ^ 3.3 objects 
(green circles) clearly indicates the superposition of fore¬ 
ground sources. The “global” escape fractions for these two 
galaxies are 140% and 187%, respectively. If we compute in¬ 
stead the “local” escape fraction only from the region with 
clear detection in the U band, we ended up with even larger 
escape fractions, of 966% and 654%, respectively. Even in 
this case the most plausible explanation is a contamination 
by a foreground galaxy prospectively close to the z ^ 3.3 
source. 

The situation is less clear for source ID=72687 in Table 
[T] This galaxy has an escape fraction of 292%, clearly not 
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Table 1. The galaxies in the VUDS-LBC/COSMOS field used to derive the LyC relative escape fraction 


IDlbc 

Field 

RAD 

(deg) 

Dec 

(deg) 

^spec 

Rmaq 

(AB) 

Umaq 

(AB) 

Misoo 

rreL 

J esc 

Cont 

13903 

COSMOS 

150.414993286 

-b2.159081936 

3.290 

22.59 

> 28.69 

-23.11 

< 0.039 

No 

74113 

COSMOS 

149.886260986 

-b2.276002169 

3.330 

23.56 

> 29.50 

-22.14 

< 0.045 

No 

56982 

COSMOS 

149.847152710 

-b2.373027802 

3.356 

23.65 

> 29.46 

-22.05 

< 0.051 

No 

72509 

COSMOS 

149.994003296 

-b2.285038948 

3.374 

23.84 

> 29.38 

-21.86 

< 0.065 

No 

16700 

COSMOS 

149.940277100 

-b2.082487583 

3.354 

24.44 

> 29.83 

-21.26 

< 0.075 

No 

61293 

COSMOS 

149.870300293 

-b2.349019766 

3.292 

23.90 

> 29.13 

-21.80 

< 0.087 

No 

21413 

COSMOS 

149.974533081 

-b2.111910582 

3.319 

24.50 

> 29.68 

-21.20 

< 0.091 

No 

31070 

COSMOS 

149.912933350 

-b2.236840248 

3.333 

24.74 

> 29.86 

-20.96 

< 0.096 

No 

64371 

COSMOS 

150.051940918 

-b2.332614422 

3.370 

24.19 

> 29.26 

-21.51 

< 0.100 

No 

49197 

COSMOS 

149.978881836 

-b2.135632515 

3.395 

24.81 

> 29.80 

-20.89 

< 0.108 

No 

58808 

COSMOS 

149.962997437 

-b2.363846302 

3.359 

24.63 

> 29.58 

-21.07 

< 0.112 

No 

68731 

COSMOS 

150.138900757 

-b2.307063818 

3.288 

24.47 

> 29.41 

-21.23 

< 0.113 

No 

32388 

COSMOS 

149.778930664 

-b2.229439735 

3.300 

24.63 

> 29.53 

-21.07 

< 0.118 

No 

16446 

COSMOS 

149.961273193 

-b2.081201315 

3.270 

24.96 

> 29.86 

-20.74 

< 0.118 

No 

69814 

COSMOS 

149.834243774 

-b2.300836325 

3.366 

24.56 

> 29.45 

-21.14 

< 0.119 

No 

3400 

Q0933 

143.354507446 

-f 28.806911469 

3.270 

24.88 

> 29.76 

-20.82 

< 0.120 

No 

1723 

COSMOS 

150.447082520 

-b2.347564697 

3.300 

23.63 

> 28.43 

-22.07 

< 0.129 

No 

56248 

COSMOS 

149.917510986 

-f2.378670692 

3.288 

24.71 

> 29.50 

-20.99 

< 0.130 

No 

63275 

COSMOS 

149.983078003 

-b2.338610888 

3.343 

24.96 

> 29.72 

-20.74 

< 0.134 

No 

72124 

COSMOS 

149.883422852 

-b2.286791325 

3.352 

24.71 

> 29.46 

-20.99 

< 0.135 

No 

51227 

COSMOS 

149.892135620 

-b2.414754868 

3.280 

24.45 

> 29.18 

-21.25 

< 0.137 

No 

50989 

COSMOS 

149.834274292 

-b2.416670561 

3.310 

24.61 

> 29.32 

-21.09 

< 0.140 

No 

51473 

COSMOS 

150.065383911 

-b2.411855936 

3.355 

25.29 

> 29.90 

-20.41 

< 0.153 

No 

8556 

Q0933 

143.382461548 

-f 28.753007889 

3.330 

24.91 

> 29.51 

-20.79 

< 0.155 

No 

17175 

Q1623 

246.469177246 

-f 26.892374039 

3.340 

24.69 

> 29.29 

-21.01 

< 0.155 

No 

47539 

COSMOS 

149.888641357 

-b2.144912004 

3.311 

25.02 

> 29.59 

-20.68 

< 0.159 

No 

30576 

COSMOS 

150.141159058 

-b2.239602089 

3.336 

25.26 

> 29.82 

-20.44 

< 0.160 

No 

10235 

COSMOS 

150.106704712 

-b2.044517517 

3.328 

25.14 

> 29.70 

-20.56 

< 0.161 

No 

70238 

COSMOS 

149.853439331 

-b2.298267603 

3.292 

24.58 

> 29.13 

-21.12 

< 0.162 

No 

42374 

COSMOS 

150.137237549 

-b2.174738884 

3.294 

25.44 

> 29.87 

-20.26 

< 0.181 

No 

19739 

COSMOS 

149.934432983 

-b2.102468014 

3.316 

25.54 

> 29.96 

-20.16 

< 0.183 

No 

12646 

Q0933 

143.328750610 

-f 28.719076157 

3.330 

25.13 

> 29.50 

-20.57 

< 0.192 

No 

10849 

Q0933 

143.360122681 

-f 28.734085083 

3.350 

25.59 

> 29.92 

-20.11 

< 0.198 

No 

65540 

COSMOS 

149.848281860 

-b2.325203896 

3.320 

25.50 

> 29.58 

-20.20 

< 0.250 

No 

4890 

COSMOS 

149.824325562 

-b2.009310722 

3.301 

25.24 

> 29.24 

-20.46 

< 0.268 

No 

19521 

COSMOS 

150.330398560 

-b2.072281122 

3.303 

22.55 

26.50 

-23.15 

0.281^“^ 

Yes 

63327 

COSMOS 

149.843582153 

-b2.340100288 

3.329 

25.02 

28.47 

-20.68 

0.446^*’^ 

No 

25155 

COSMOS 

149.857284546 

-b2.270117283 

3.364 

25.24 

28.61 

-20.46 

0.479(‘^^ 

No 

40998 

COSMOS 

149.952316284 

-b2.182848215 

3.321 

25.19 

28.36 

-20.51 

0.577^'*) 

Yes 

48734 

COSMOS 

149.815368652 

-b2.138106108 

3.384 

24.99 

27.99 

-20.71 

0.674^'^) 

Yes 

47357 

COSMOS 

150.096496582 

-b2.145905495 

3.292 

24.91 

27.78 

-20.79 

0.760(-''^ 

Yes 

52336 

COSMOS 

149.980270386 

-b2.405712605 

3.371 

24.85 

27.06 

-20.85 

1.400 

Yes 

55594 

COSMOS 

149.939086914 

-b2.388362885 

3.334 

24.17 

26.06 

-21.53 

1.873 

Yes 

16669 

COSMOS 

149.994964600 

-f2.080330372 

3.367 

24.39 

26.06 

-21.31 

2.295 

Yes 

72687 

COSMOS 

149.877334595 

-b2.283946991 

3.328 

24.14 

25.55 

-21.56 

2.920 

Yes 


The galaxies have been sorted in ascending order of relative escape fraction. The upper part of the table summarizes the non 
detections (i.e. S/N < 1) in the LBC U-band, for which an upper limit to the relative escape fraction has been provided. The 
central sub-sample includes galaxies with apparent Lyman Continuum detection, while the last four galaxies have escape fraction 
larger than 1.0 (i.e. 100%) and are thus considered as unreliable and probably due to foreground contamination, as discussed in 
the main text. The last column, labeled Cont, indicates whether a galaxy is possibly contaminated by foreground emission. Notes 
on individual objects: (a) local escape fraction 1.75, calculated on the LyC emission region only; (b,c) they can be plausible LyC 
emitters; (d) local escape fraction 4.51; (e) local escape fraction 6.40; (f) local escape fraction 5.13. 


physical under reasonable assumptions for the stellar popu¬ 
lations involved and for the IGM absorption. From an anal¬ 
ysis of the HST image (FigHI we cannot confirm with high 
confidence the presence of two different objects as in the 
previous cases, although two distinct blobs are seen in the 
HST thumbnail. From a detailed analysis of the LBC im¬ 
ages, we find a small shift of ~ 0.5 pixel, corresponding to 
~0.1 arcsec, from the center of the source in the R band 
and in the U band. The separation between the two blobs in 
the HST image is slightly larger, of 0.3 arcsec. This could 


be an other case of blending of a z 3.3 galaxy with a 
foreground source, but based on the present data we can¬ 
not exclude that the two blobs detected by HST are all at 
z > 3, and this could be an example of a stellar-powered 
LyC emission. A detail analysis of the VUDS spectrum in¬ 
dicates that the redshift of this source is based only on a 
noisy Lyman-a and strong CIV and SilV lines, all in ab¬ 
sorption. There are also other absorption lines not associ¬ 
ated with the z = 3.32 solution, but it was not possible 
to derive a secure redshift from these lines. Moreover, this 
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Fig. 1. The cutouts of source ID=16669 in the U band (left) and in the R band (center) by the LBC instrument. The 
cutout on the right shows the same sky area observed by HST in the I814W filter. The position of the galaxy with 
spectroscopic redshift is indicated by the green circle. A diffuse emission at the north of the z ^ 3.3 galaxy is detected 
in all the three bands. The red ellipse in the central panel indicates the area covered by source ID=16669 in the LBC R 
band. This region has been used to compute the “global” escape fraction provided in the last column of Table [TJ The 
size of each image is 20 arcsec by 20 arcsec. 


source has strange colors (u — g = 0.13, g — K = 2.87 
and K — 8.0p,m = 1.94), which are not typical of a z ~ 3 
st ar forming galaxy. I nde ed, the Ultravista DR l catalogs 
by Ilbert et al. (2013) and Muzzin et al. (2013) give a pho- 
toinetric redshift of 0.596 and 0.628, respectively, for this 
object. Further analysis of this very peculiar source will be 
crucial to understand its nature. 


4.1.2. Galaxies with < 100% but detected LyC emission 

In Table [I] there are 6 additional galaxies with detection in 
the LyC wavelengths and estimated values of < 100%. 

A careful analysis of the LBC images in the U and 
the R bands, complemented by the HST imaging in the 
I814W band, allows us to distinguish between real LyC 
emission and foreground contamination. Four (ID=19521, 
40998, 48734 and 47357) out of these six sources show a sit¬ 
uation similar to the four previously analysed galaxies: the 
high spatial resolution images by HST ('Fig lA.2|) show the 
presence of possible faint contaminants with spatial offsets 
from the spectroscopically identified galaxies at z ~ 3.3. 
To check if they are foreground galaxies, we adopt, as in 
the previous cases, the test of the “local” escape fraction 
to verify if it is consistent with the emission from plausi¬ 
ble stellar populations and with the IGM attenuation. The 
escape fraction has been computed restricting the area to 
the region covered by the LyC emission in the U band, in¬ 
stead of computing the fluxes in the region defined by the 
non-ionizing band (in this case the LBC R-band). For all of 
these four galaxies, the escape fraction computed “locally” 
on the U band emission further exceeds the physical limit 
of 100%, in some cases reaching 600% or more (see notes 
to Table [T|). 

For the remaining two galaxies (ID=25155 and 63327) 
the presence of a contaminating foreground is not obvi¬ 
ous from HST data, and their “local” escape fractions are 
< 100% if we limit the analysis to the region with U band 
emission. The galaxy 25155 (Figl3]) is isolated and only 
checking carefully the HST I814W image there is a hint of 


a very faint blob slightly offset from the z = 3.3 galaxy. It 
has a distance of 0.1 arcsec from the bright object, thus 
it is very hard to draw any conclusions on its true na¬ 
ture. We will assume in the following that it could be a 
genuine LyC emitter at z ~ 3.3. This galaxy has a rela¬ 
tively faint magnitude U = 28.6 and it is marginally de¬ 
tected in the LBC COSMOS U band image at S/N = 2.2, 
with a resulting escape fraction of ^ 48%. Its optical spec¬ 
trum by the VUDS survey shows typical absorption lines 
of high-z galaxies (Ly-a, SilV, CIV) so it could resem ble 
th e “lonl” source discus sed in Vanzella et al. (2012a) and 


Vanzella et al. (2015) 


In addition, also the detected U-band flux for galaxy 
ID=63327 cannot be clearly explained by the presence of 
a lower redshift object contaminating the z ^ 3.3 source’s 
line of sight. We have thus another example of a plausi¬ 
ble ionizing source by stellar radiation, with = 44.6% 
and a marginal detection at 17 = 28.5 at S/N = 2.2. FigH] 
shows this source in the LBC U band (left), in the LBC R 
band (center) and in the I814W band by HST (right). The 
HST image shows two blobs separated by ^0.7 arcsec, but 
even at this high resolution we cannot judge whether this 
is a unique object or two distinct galaxies. There is still the 
possibility that there is a foreground contaminant roughly 
aligned with the z ~ 3.3 galaxy making an apparent blend 
that is indistinguishable even from HST data. For further 
exploring this possibility we used the high resolution of the 
HST I814W image to carry out a de-blended photometry 
both in the U and in the R band with ConvPhot, as de¬ 
scribed in the previous sections, but now adopting the HST 
image as a prior. We obtain an upper limit of 17 > 29.4 at Itr 
for the northern blob, while for the southern region we get 
a marginal detection U = 28.9 at 2cr. Translating the latter 
value into an escape fraction, we obtain f/fl = 119%, thus 
indicating that this can be a lower redshift interloper, sim¬ 
ilar to the ones discussed in the previous section. However, 
since this value of the escape fraction can be lower than 
100 % for slightly bluer stellar populations or a more trans¬ 
parent line of sight, we still consider this object as a candi- 
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Fig. 2. The cutouts of source ID=72687 in the U band (left) and in the R band (central panel) by the LBC instrument. 
The right cutout shows the same sky area observed by HST in the I814W filter. The position of the spectroscopic redshift 
is indicated by the green circle. The size of each image is 14 arcsec by 14 arcsec. In this case the presence of two distinct 
objects is not secure, despite the availability of the HST image, but there is some hints of the presence of two blended 
galaxies in the latter. The object in the U band is slightly offset (~ 0.1 arcsec) from the center of the source detected in 
the LBC R and HST 1814W bands. 


date LyC emitter, requiring further analysis in the future. 
We have checked also the VUDS spectrum of this source 
looking for possible absorption features not consistent with 
the redshift of the z ~ 3.3 object, but the S/N is too low 
to carry out this test. Its optical spectrum shows only a 
relatively narrow Lyman-a in emission, with a symmet¬ 
ric profile. Deep UV images with HST or deep NIR spec- 
tro scopy could finally reveal its nature, as recently shown 
by Siana et al. (2015) for similar cases in the SSA22 area. 

Summarizing, we have found evidence of contami¬ 
nation in 8 out of 45 cases (18%). This fraction is 
slightly larger than, but still consistent with the ex- 
pec ted fraction of contami nation (12-15%) derived follow¬ 
ing Vanzella et al. (2010a) assuming a depth of 17 = 28.5 — 
29.5 and a seeing of 1.0 arcsec. We can thus conclude that 
there is a high probability that these 8 galaxies detected in 
the U band are all foreground interlopers. 


5. Results 


5.1. The escape fraction derived on stacked images 

We derived the stacked images in the U and R bands for 
the 37 galaxies presented in Table [T] marked with the la¬ 
bel Cont = No, excluding the 8 galaxies where we have 
an evidence of possible foreground contamination, in order 
to measure the mean escape fraction of the population of 
star forming galaxies. We have included in the stack the 
two galaxies (1D=25155 and 1D=63327) with marginal de¬ 
tection in the U band. We adop ted the same technique de¬ 
scribed in Boutsia et al. (2011) In summary, we extract a 
thumbnail of ^ Abarcsec x 45arcsec centered around the 
0 ~ 3.3 galaxies in the U and R band images by LBC, de¬ 
tect all the sources in the R band and mask all of them in 
both bands, except for the central target. Then, we aver¬ 
age all the thumbnails of the 37 galaxies using the inverse 
square RMS map as a weight. No clipping has been ap¬ 
plied during the stacking. The resulting images (zoomed) 
are shown in Figl^ It is immediately clear the high S/N 
in the stacked R band image, given that each individual 
source is detected at S/N > 10 in the original frame. It is 



Fig. 5. The stack of the 37 galaxies with 3.27 < Zspec < 
3.40 in the U (left) and R bands (right) for the VUDS- 
LBC/COSMOS, Q0933, and Q1623 fields. Other objects 
detected in the R band but not associated with the galaxies 
at z 3.3 (at the center of the FITS thumbnail) have been 
masked both in the U and in the R band image using the 
segmentation map derived from the R band. The size of 
each image is 14 arcsec by 14 arcsec. 


also evident that no detection is apparent in the U band 
stack. 

To be more quantitative, we adopted the same method 
described in the previous section for the individual sources, 
namely we used the ConvPhot software to carry out a care¬ 
ful evaluation of the depth of the stacked images. The 
stacked source in the R band has a mean magnitude of 
R = 24.56 with a S/N = 283. On the other side, for the U 
band we can derive only an upper limit to the emitted flux, 
reaching U = 31.4 at S/N = 1 (or equivalently U = 30.2 
at S/N = 3). This translates into a lower limit for the 
flux ratio of (Fn/Fu)obs = 545.1, corresponding to an up¬ 
per limit to the relative escape fraction of //// < 2.0% at 
Icr, or equivalently 6.0% at the 99.7% confidence level. If 
we exclude from the stack the two sources 1D=25155 and 
1D=63327, i.e. the galaxies detected individually in the U 
band at ^ 2a level, the previous results are not affected 
at an appreciable level. The stack of the 35 sources indeed 
gives a limit to the flux ratio of 549.5 corresponding to a 
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Fig. 3. The cutouts of source ID=25155 in the U band (left) and in the R band (center) by the LBC instrument. The 
cutout on the right shows the same sky area observed by HST in the I814W filter. The position of the spectroscopic 
redshift is indicated by the green circle. The size of the cutouts is 4.3 arcsec. This object seems isolated in the HST 
I814W band, but there is a possible hint for the presence of a very faint source on the north-west of the bright galaxy at 
z ^ 3.3. There is a detection at ^ 2a level in the U band of LBC. 



Fig. 4. The cutouts of source ID=63327 in the U band (left) and in the R band (center) by the LBC instrument. The 
cutout on the right shows the same sky area observed by HST in the I814W filter. The position of the spectroscopic 
redshift is indicated by the green circle. The size of the cutouts is 15 arcsec. In the HST image there is a faint blob 
towards the south-west, but it is blended with the galaxy at z ~ 3.33 in the LBC R band. In the LBC U band there is a 
very faint emission slightly offset towards the south w.r.t. the R band position. 


similar limit for the escape fraction of 2.0% at Icr. Thus 
the inclusion of two faint galaxies with marginal detection 
in the LyC does not modify the actual limit reached by the 
stacking. 

We split the whole sample in different sub-samples. We 
explore the escape fraction for the bright and faint sub¬ 
samples. Considering the 11 galaxies with R < 24.5 within 
the total sample of 37 objects used above, we end up with 
a weighted-mean magnitude of R = 23.93 and an upper 
limit for the escape fraction of 2.4% at Icr. The stack of 
the remaining 26 sources with R > 24.5 gives instead a 
weighted mean of i? = 24.92 and < 3.1% at Icr (or 
equivalently < 9.4% at 3 ct level). In this case the 
slightly worse limit to the LyC escape fraction is only due 
to the faintness of the sub-sample. 

Figin] shows the escape fraction at 1 tr as a function of 
the observed R band magnitudes for all individual sources 
and for the stacks. The trend of larger with fainter R 
band magnitudes is simply due to the depth of the LBC 
image in the U band (U=29.7 AB mag at S/N=l), and it 


Table 2. The escape fraction of stacks 


Sample 

Ngal 

Rmag 

(AB) 

nrel 

Jesc 

(la) 

bright {R < 24.5) 

11 

23.93 

<0.024 

faint {R > 24.5) 

26 

24.92 

<0.031 

all sources 

37 

24.56 

<0.020 


is not due to a physical trend in the galaxy population (i.e. 
higher escape fraction for fainter sources). 


The results of the stacks are summarized in Tableland 
are also plotted in Fig|B]with magenta, red, and blue arrows 
for the total sample and the bright and faint sub-samples, 
respectively. These limits are a factor of 2.5 more s trin- 
gent than the previous results by Boutsia et al. (2011) and 
demonstrate the power of stacking a large sample of galax- 
ies to reach very deep constraints to the ionizing escape 
fraction at z ^ 3.3. 
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Fig. 6. The measured values and the upper limits at 1 ct for the LyC relative escape fraction of galaxies at z ^ 3.3 
in the COSMOS field. Blac k descending arrows sh ow the upper limits associated to individual galaxies, while the cyan 
arrow is the limit derived in Boutsia et al. (2011) by stacking 11 galaxies. The red arrow is the limit derived using only 
relatively bright galaxies with R < 24.5, corresponding roughly to L > 1.2L* (M 1500 < —21.2). The blue arrow shows 
the upper limit for fainter galaxies (stack of 26 objects), while the magenta one indicates the upper limit for the whole 
sample (stack of 37 galaxies). The magenta dashed line is the expected for a limit in the U band of 29.7 magnitudes 
(at S/N=l, assuming an intrinsic ratio L 1500 /L 900 = 3 and a mean IGM transmission of 0.28). This shows that the 
trend of larger escape fractions with fainter R-band magnitudes is due to the depth of the LBC image in the U band, 
and it is probably not due to a physical trend in the population of star-forming galaxies. Asterisks show the galaxies 
with detection in the U band and a significant LyC emission. The grey hexagons and the upper arrows indicate likely 
foreground contamination, where the global or the local escape fractions exceed 100%. Only two sources (ID=25155 and 
ID=63327, black asterisks without the grey hexagons) are an example of possible LyC emitters. 


5.2. UV Background from star-forming galaxies at z 3 

Using the results obtained in the previous section, we can 
now derive a constraint for the UV background produced 
by the population of star-forming galaxies at z ^ 3. We 
derived the luminosity density at 900 A rest frame as in 
Boutsia et al. (2011)[ starting from the analogous quantity 
at 1500 A rest and correcting it for the intrinsic slope of 
the spectrum of these galaxies (adopting L 1500 /L 900 = 3) 
and the derived limit for the escape fraction of 2%: 

P900 = P1500 • fefc ■ {L90o/Li50o)int ■ (3) 

As it can be derived from equation 1, the product 
fesciLgoo/ Li5oo)int IS equivalent to (Fu/Fii)obsexp{TiGM) 


which are measured quantities less subject to systematic 
uncertainties or assumptions. 


The only parameter affecting these calculations is the 
input luminosity density of non-ionizing radiation, pisoo- 
To obtain this value it is sufhcient to integrate an observed 
LF, multiplied by the absolute luminosity, down to a spe¬ 
cific m agnitude limit. We decided to adopt the parametriza- 
tion of Reddy & Steidel (2009) for the LF at z ^ 3 after 
converting it fr om 1700 to 1500 A re st frame wavelengths, 
as described in Boutsia et al. (2011) We adopt a slope of 
a = —1.73 and integrate down to an absolute magnitude of 
Miboo = —20.2, which is equivalent to L ~ 0.5L*(z = 3). 
We choose this limit since it corresponds to the limiting 
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magnitude of our sample, R 


nosity density is pisoo = 1-27 x 10^®erg s ^ Hz ^ Mpc ® 


25.5. The resulting lumi- 
erg s~ 

Adopting an upper limit for the escape fraction of < 2.0%, 
we obtain a limit to the ionizing emissivity of pg^Q < 
0.85 X lO^^erg Hz~^ Mpc~^ at Icr level, or equivalently 
Pqoo — 2-55 X lO^'^erg Hz~^ Mpc~^ at 3cr. 

We can translate this ionizing luminosity density into 
a meta-galactic hydrogen photo-ionization rate r_i 2 , (ion¬ 
ization rate of HI in unit of 10~^^s~^atom~^) using the 


formula by Meiksin (2009) and Mostardi et al. (2013) 


-12 — 


10^" • Ploo ■ ■Al-{1 + zf 

hp ■ {3 + \auv\) 


(4) 


where api = 6.3 x 10“^®cm^ is the hydrogen photo¬ 
ionization cross section at A = 912 A, apy = —3.0 
is the spectral slope of th e stellar ionizing radiation 
(from Mostardi et al. (2013)1 and hp is the Planck con¬ 
stant. T'he mean free path of ionizing photons, Al is 83.5 
proper Mpc at z=3.3, according to recent estimates by 
Worseck et al. (2014a) The Icr uncertainty on Al at z = 
3.3 is 4.8 Mpc (i.e. 6%). The uncertainties on the emissivity 
of galaxies at 1500 A rest frame (due to the uncertainties 
on the parametrization of the galaxy LFs) are thus domi¬ 
nating the source of errors on r_i 2 , as shown in Table 3. 
We will discuss in the next section the implications of dif¬ 
ferent parametrization of the galaxy LFs on the HI photo¬ 
ionization rata _ 

Following Faucher-Giguere et al. (2009)[ at z ~ 3 the 
galaxy radiation can be considered rather uniform, and in 
this case the net effect of discrete absorbers is a hard¬ 
ening of the emissivity, due to the IGM filtering. As a 
consequence, a correct val ue for galaxies can be apy = 
— 1.8 (from Appendix D of Faucher-Giguere et al. (2009)1, 
which is in agreement with the value provided By 
Haardt & Madau (2012)|(see their Fig. 10). 

We thus obtain an upper limit F_i 2 < 0.12 at Icr 
(or 0.36 at 3cr) for the ionization rate of relatively 
bright (Mi 5 oo < —20.2, or L > 0.5L*) star forming 
galaxies at the mean redshift of our sample (z=3.3). For 
comparison, a recent estimate of the ionization rate by 


Becker & Bolton (2013)| gives r_i 2 = 0.79l°'^® at z=3.2, 
Figlll 


well above our limit. FiglTjshows the derived upper limit (at 
99.7% c.l. or 3cr) for r_i 2 from the VUDS-LBC/COSMOS 
data, compared with recent determinations of the UV 
background (UVB) in the literature. These estimates 


properties (JBolton et al. (2005) 

IKirkman et al. (2005) 

Faucher-Giguere et al. (2008) 

Wyithe & Bolton (2011) 

Becker & Bolton (2013)) or by the proximity effect 

(Calverley et al. (2011)|). The data shown in FiglTJ derived 

from hux decrement analysis in 

the IGM Lyman forest 


of QSO spectra (Bolton et al. 2005; Faucher-Giguere et 
al. 2008; Wyithe & Bolton 2011), have been scaled to the 
same cosmology and to the same IGM temperature-density 
relation, as described in Calverley et al. (2011) and in 
Giallongo et al. (2012). It is clear from this plot that 
the ionizing emission of relatively bright (L > 0.5L*) 
star forming galaxies is insufficient to keep the Universe 
reionized at z ^ 3.3 at more than 99.7% probability. 
It is worth mentioning that the recent determinations 


of Becker & Bolton (2013) take into account a series of 
systematic effects on the derivation of r_i 2 , and thus their 
error bars can be considered as a robust estimates for the 


0.5 


7 - 0.5 - 


-1 


- 1.5 



■ Bolton et al. (2005) 

* Kirkman et al. (2005) 

• Faueher-Giguere et al. (2008) 

•Wyithe and Bolton (2011) 

• Calverly et al. (2011) 

* Becker and Bolton (2013) 

^PVUDS/LBC L>0.5L* Galaxies (3cj) 


4 

Redshift 


Fig. 7. The evolution of the UVB, (ogio(r_ 72 ), its a func¬ 
tion of redshift. These measurements are typically derived 
by an analysis of the Lyman forest properties or by the 
proximity effect. The 3cr upper limit from this work (r_i 2 < 
0.36) is indicated by the magenta arrow. This limit, ob¬ 
tained by star forming galaxies with L > 0.5L* at z ^ 3.3 
shows that the bright galaxies alone cannot be responsible 
for the UVB at this redshift. 


uncertainties on the UVB. Our upper limit is clearly lower 
than their minimum envelope, giving convincing argument 
against a major contribution of bright star forming galaxies 
to the ionizing emissivity at z=3. In the next section we 
will discuss in detail the possible way outs for the missing 
HI ionizing photons at z ~ 3. 


6. Discussion 

Our results are based on two main assumptions, namely 
the mean IGM transmission < ea;p(— t/gm) > and, in a 
more subtle way, the adopted value for the intrinsic ratio 
(Li 50 o/L 9 oo)mt- We will discuss in the following sections 
the implications of the choice of these two parameters on 
our conclusions. 


6.1. The variance of the IGM and ISM transmissions 

In the previous sections we have adopted a mean value of 
0.28 for the mean IGM transmission close to the Lyman 
limit. This value however is subject to large cosmic vari¬ 
ance due to differences on the absorption patterns in vari¬ 
ous lines of sight, especially due to the patchy nature of the 
Lyman limit systems (LLS) and Damped Lyman-a systems 
(DLAs), which are able to block the path of ionizing pho¬ 
tons through the IGM. The IGM absorption at 900 A is 
thus very stochastic because it is related to the probability 
of intercepting dense absorbers, with log{NHi) > 17. 

We have estimated the variance of the IGM transmis¬ 
sion exp{—TicM) at ^ 900 A rest frame adopting the Monte 
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Carlo simulations of Inoue&Iwata_^200^ and th e recent 
up dated absorber statistic s by Inoue et al. (2014) used also 
by Vanzella et al. (2015) Recently, Thornas et ai. (2015) 
estimated the same quantity using deep VUUS spectra 
of galaxies in the redshift interval 2.5 < z < 5.5, find¬ 
ing an excellent agreement with the IGM transmission ob¬ 
tained from quasars’ lines of sight at z < 4. We consid¬ 
ered 10,000 different realizations of the IGM transmission 
at z ~ 3.3 and convolved each of them with the profile 
of the LEG U band filter, thus obtaining a mean trans¬ 
mission and its variance. We obtained a mean value of 
< exp[—TiGM) > = 0.25, similar to the a dopted value of 
0.28 based on the Worseck et al. (2014a) IGM transmis¬ 
sion, and a dispersion of ^ 0.2. The resultirig distribution of 
IGM transmission at z ^ 3.3 is asymmetric, with a tail to¬ 
ward low exp{—TiGM) due to the presence of DLAs and LL 
systems. A fraction of 11% of all the simulated lines of sight 
shows a transmission exp{—TiGM) less than 0.01, thus it is 
important to have a large number of galaxies with measured 
LyC escape fraction in order to minimize the stochastic ef¬ 
fect of intervening absorbers. Considering our case of 37 
different lines of sight, the dispersion of < ea;p(— t/gm) > 
is 0.03. 

Another source for the skewed probability distribution 
function towards low values of the escape fraction of galax¬ 
ies is the patchiness of the ISM, which has been recently ex¬ 
plored by Cen & Kimm (2015) Using high-resolution large 
scale simulations with radiative transfer calculations, they 
show that LyC photons can escape from a galaxy through 
small solid angles with complete transparency, while the 
majority of the remaining LoS are opaque. Adopting their 
simulations at z = 4 and at a virial mass 10^° < M < 
10^^Mq, we estimate that for an input escape fraction of 
2% with 37 galaxies the confidence interval at 1 cr level is 
between = 1.4% and 2.7%. Simulations are not avail¬ 
able for our redshift and virial mass range, but the trend 
with redshift and masses derived by Cen & Kimm (2015) 
shows that our extrapolation at z 3.3 and higher lu- 
minosities should be fair enough. Combining together the 
uncertainties due both to the stochasticity of the IGM and 
the variance of the ISM in different LoS, the confidence 
interval at 1 cr is between = 1.3% and 2.8%. 

Thus we can conclude that with our relatively large sam¬ 
ple of galaxies we are able to limit the variance due to 
both IGM transmission and patchy ISM of the star form¬ 
ing galaxies. We can exclude that the lack of ionizing pho¬ 
tons we observed in the VUDS-LBC/COSMOS sample is 
due to the low number statistics and variance in the IGM 
or ISM transmission. In the near future, the combination 
of HST imaging, deep UV observations and large spectro¬ 
scopic surveys will provide more than 100 galaxies where 
these measurements can be repeated. 

6.2. The intrinsic ratio (Lisoo/Agoo)™* 

Our choice of the (Li 5 oo/A 9 oo)mi ratio does not enter 
directly in the evaluation of the total ionizing emissiv- 
ity since the dependencies of Pqqq from the intrinsic ratio 
{Li5oo/L goo)int is Canceled out, once equations 1 and 3 are 
combined. As we have already discussed, in fact, the value 
of the HI photo-ionization rate does depend only on the 
observed ratio between the U and the R band photometry 
and on the IGM correction. However, the (Lisoo/Agoo)™* 
ratio can be used to provide us a with a physical interpre¬ 


tation of the amount of UV ionizing photons escaping from 
a galaxy, and enters in our definition of the sample. We will 
discuss this latter aspect here. 

The relative escape fraction derived for our sample of 
star forming galaxies assumes a given ratio between the in¬ 
trinsic luminosity at 1500 and 900 A rest frame, here taken 
to be (Li 5 oo/A 9 oo)mt = 3. This value has been adopted for 
a direct comparison with previous works dealing with sim- 


ilar subjects (Steidel et al. (2001) Shapley et al. (2006)). 


The intrinsic ratio (Li 5 oo/Lgoo)mt depends on the star for- 
mation history, stellar populations, age, metallicity, initial 
mass function of the investigated galaxies. 


(L 


Vanzella et al. (2010b) showed that the intrinsic value 
1500 /Lgoo) ini Can be as large as 7 if considering mean 
values for the physical parameters of high-z star-forming 
galaxies. Adopting this new value would approximately 
double the escape fraction we measured for the VUDS- 
LBC/COSMOS sample, but would not change our con¬ 
clusions on their contribution to the ionizing emissivity. 
Moreover, adopting a value of 7 for (Li 5 oo/A 9 oo)int would 
result in even higher values for the relative escape frac¬ 
tion of the galaxies which already show > 100%, thus 
strengthening our conclusions about their nature of low red¬ 
shift interlopers. 


Alternatively, Topping & Shull (2015) adopted a recent 
rendition of stellar evolutionary tracks as input in their 
galaxy synthesis model, finding a 50% increase in the 
LyC production efficiencies w.r.t. previous calculations. 
Changing the intrinsic value of (Li 5 oo/Lgoo)mt to 2 would 
not change our conclusions on the ionizing emissivity, as dis¬ 
cussed above. Moreover, since the relative escape fraction 
we derived with the “local” method is significantly greater 
than 100% (175-966%) for all the galaxies that were con¬ 
sidered to be affected by foreground contamination, a cor¬ 
rection of 50% downward would only partially alleviate the 
tension for one object (ID=19521), but would not solve the 
problem for the other galaxies. 

We can thus conclude that the assumed value of 
(Ai 5 oo/A 9 oo)mt = 3 is not significantly affecting our con¬ 
clusions on the ionizing emissivity of star forming galaxies 
at z ~ 3.3. 


6.3. Comparison with other results 


According to Cooke et al. (2014), the LBG selection could 
be biased against LyC sources, since a galaxy with signifi¬ 
cant emission at A < 912 A rest-frame would avoid the clas¬ 
sical drop-out selection. We show here that our sample is 
not affected by these systematics. In fact, even assuming an 
escape fraction of 100% and a completely transparent IGM 
(very unlikely), the typical color of a galaxy at z ^ 3 would 
he U—R > 1.2, above the color criteria normally adopted to 
select ca ndidate galaxies for spe ctros copic follow up (see fo r 
example Le Fevre et al. ( 2015)| and [Boutsia et al. (2014) ). 
A recent paper (jVanzella et al. (2015)P confirms our con- 
clusions with st ate-of-the-art simula tions adopting the up¬ 
dated recipes of Inoue et al. (2014) for the IGM transmis¬ 
sion. 

The result found from the VUDS-LBC/COSMOS 
sample can be compared with recent achievements on the 


detection of LyC emitters at z > 3. 

Iwata et al. (2009) 

Nestor et al. (2011)1 

Nestor et a 

L (2013) and 

Mostardi et al. (2013) found a significant number of 


ionizing sources m their surveys. In the 


area, 
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Iwata et al. (2009) found evidence for /JJ)! > 4% for the 
whole population of 73 LBGs analysed, but provided 
also evidence for a significant spatial offset between the 
ionizing emission and the non-ionizing one for some of 
them. Moreover, there are few objects detected in LyC 
with a flux ratio L’lsoo /T^goo smaller than the typical values 
expected from PopII stellar population with standard 
Salpeter IMF (or other non top-heavy ones), even assum¬ 
ing a transparent IGM. Moreover, they didn’t carry out a 
local estimate of the escape fraction in the position of the 
detected LyC sources. 

On the same SSA22 
show 


field, 


Nestor et al. (2011) and 
Nestor et al. (2013)| show nine LBGs and 20 Lyman-a 
emitters (LAEs) with LyC detection out of a sample of 
41 LBGs and 91 LAEs (all spectroscopically confirmed). 
They started from a different narrow band image cen- 
ter ed at ^ 3640 A w hich is deeper than the one used 
by Iwata et al. (2009)| at ^ 3590 A. A careful analy¬ 
sis of their LBC detections, however, shows that even in 
this case the LyC emission for many z ~ 3 sources is 
offset by 0.4-1.0 arcsec. The observed ratio between the 
900 and 1500 A rest frame emission is difficult to recon- 
cile with what expecte d by standard stellar populations 
(Vanzella et al. (2012a)). The HST images in I814W filter 
available for few of them show the presence of clearly sep¬ 
arated galaxies, sometimes fainter at 1500 A than in the 
ionizing continuum (i.e. their C16 object), with a resulting 
escape fraction well exceeding 1000 % if estimated in the 
LyC position. For the majority of them no HST imaging 
is available but even ground-based images often show the 
presence of slightly offset emission in LyC w.r.t. the non ion- 
izing continuum. Simil ar conclusions can be reached for the 
Mostardi et al^_(^0jj} sample where they adopt the same 
analysis as in |Nestor et al. (2013)| They found at z ~ 2.8 
four LyC emitters out of 49 LBC galaxies and 7 LyC emit¬ 
ters out of 91 Lyman-a emitters. In this case the lack of 
high spatial resolution data from HST for the majority 
of the sample prevents any detailed analysis about pos¬ 
sible contamination by interlopers/foregrounds. These con- 
clusions have been strengthened by recent observations by 
Siana et al. (2015)| where they found no convincing detec- 
ir deep H£ 


tion in their deep HST imaging with UVIS of 5 LyC emit- 
ters extracted from the sample of Nestor et al. (2011)[ or by 


Mostardi et al. (2015)[ where only one robust LyC emitter 


ha s been found after a re-a nalysis of a sample of 16 galaxies 
by Mostardi et al. (2013) Mor e interestingly, the only ro- 
bust candidate LyC emitter by|Mostardi et al. (2015)| the 
galaxy MD5b, has an observed ratio ^uv/^LyC = 4.0±2.0, 
equivalent to a relative escape fraction of 75%, assuming 
complete transmission of the IGM. If instead a mean value 
of < e xpj—TiGM) >= 0-4 at z ~ 3.1 will be adopted, fol¬ 


lowing Inoue et al. (2014) the relative escape fraction of 
MD5b turns out to be 188%. Imposing the constraint of 
a physical value for the relative escape fraction of < 
100%, the line of sight of MD5b must be very transparent 
exp[—TiGM) > 0.75, which corresponds to a probability 
< 10“'*, according to Inoue et al. (2014) This could be an 


In summary, it is difficult to compare in a homogeneous 
way different results obtained through imaging of differ¬ 
ent quality in terms of depth and angular resolution. It 
appears however that the apparent discrepancies from dif¬ 
ferent databases tend to converge toward low values for 
the average escape fraction if we remove from the analy¬ 
sis the probable low-z contaminants, i.e. the objects with 
flux ratios F 900 /F 1500 significantly higher than that pre¬ 
dicted by spectral synthesis models, where a spatial offset 


between Fqqq and f isop is present ( [Vanzella et al. (2012a) 
Siana et al. (2015) I. 


6.4. The sources of ionizing photons at z > i 

In the previous section we have shown that bright 
(L > 0.5L*) galaxies at z ^ 3.3 cannot provide alone 
the required LyC emissivity (F_i 2 < 0.36 at 99.7% 
confidence level) to explain the measured UVB. The ACN 
contribution at z = 3.3 is r_i 2 = 0.37 (from Fig. 8 a of 
Haardt & Madau 2012). The measured value at z = 3.2 is 
r_i 2 = 0 . 79 to;i 9 from Becker & Bolton (2013). Even con¬ 
sidering the uncertainties, there is a tension between the 
measured value of r_i 2 and the summed contributions of 
ACNs plus bright galaxies at z = 3.3. This implies that an 
additional source of ionizing photons (r_i 2 = 0.06 — 0.33) 
is probably required in order to reach the observed value of 
the photo-ionizing background, taking into account its un¬ 
certainties. Looking for ionizing sources among the fainter 
galaxy population is the obvious next step. The faint pop- 


escape fraction I 

Razoumov & Sommer-Larsen (2010) 

Alvarez, Finlator & 1 

.'renti (2012) Fontanot et al. (2014) 

Faardekooper et al. (2013)1 

Kimm Cen (2014) 

Wise et al. (2014) Faardekooper et al. (2015) 

Roy et al. (2015) IDuncan V Conselice 201511. In the 


the escape fraction, slope a of the LF and magnitude 
limit Mi 5 oo, adopted to compute the ionizing emissivity, 
looking for the sources responsible for the observed UVB 
at z ^ 3.3. A summary is provided in Table (3] 

Adopting quite steep LFs as derived by 


Reddy & Steidel (2009) with a power-law index a = —1.73, 
and integrating down to M 1500 = —17.2 (or R = 28.5), we 
obtain an ionization rate of F_i 2 < 0.38 at 1 cr, assuming a 
constant escape fraction of 2 %, as for the brighter galaxies. 

The relative escape fraction of galaxies brighter 
than Misoo —20.2 (L = 0.5L*) should be ^ 16% in 
order to provide a photo-ionization rate r_i 2 ^ 0.95, 
in agreement with the observed value at z 3 shown 
in Figl3 Alternatively, the escape fraction can remain 
of the order of few percent but the LF must be inte- 
grated down to very fa int magnitudes, at M isqq = —14.2 

have 


(Fontanot et al. (2014)[ 


1^2 

servative. 


0.9. The 
since we are 
.73). At z ~ 


Alavi et al. (2014) I 
calculations 


above are 
assuming a rather 

the estimated slopes of the 


to 

quite con- 
steep LF 


indication that also this galaxy is a low-z contaminant, sim- galaxv LF are generally flatter, with a ~ 1—1.4.—1.5 

ilar to the other cases studied by Mostardi et al. (2015) (Sawiki & Thompson (2006)1 

Finkelstein et al. (2014) 

If we compare the apparent statistic of LyC detections 

Farsa et al. (2015)P, with the exception of the Lh' o 


ies (22%), similar to the value by Nestor et al. (2011)] 
significantly more than the value by Mostardi et ai. (21 


but 


013) 


( 8 %), probably due to their shallow narrow band imaging. 


— 1.8 at z > 3. Adopting the parametrization of 


a — i.» at z _ _. _ 

Sawiki & Thompson (2006) | with a = —1.43, and inte- 
grating down to AZisoo = —14.2, the ionizing emissivity is 
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r -12 £ 0.31, insufficient to ke ep the Universe reion ized. 
Similarly, if we adopt the LF by ' 
a slightly steeper slope {a = —. 


Cucciati et al. (2012) with 

~my, 


but brighter M*, we 
find an emissivity of r_i 2 < 0.26 at M 1500 = —14.2, far 
from providing enough ionizing photons to maintain the 
Universe reionized. The latter LF is based on the VVDS 


(Le Fevre et al. (2013)1 spectroscopic survey, so it should 
be robust against int erlopers in the z ^ 3 LF. At z > 3, 
Cassata et al. (2011) find that the LF of LAEs is relatively 
steep, with a slope of ~ —1.8. However, these galaxies 
are much fainter than the typical UV selected galaxies 
discovered at the same redshifts. They measured a star 
formation r ate density which is a f actor of 5 lower than 
the one by |Reddy & Steidel (2009)1 with the implication 
of a lesser contribution of the LAEs to the HI ionizing 
background at z ~ 3. Table [3] summarizes the different HI 
photo-ionization rates derived for different parametrization 
of the z ^ 3 galaxy LF and at different luminosities. 

As a _ comparison with other works 

(e-g- 


_ Iwata et al. (2009) 

Mostardi et al. 


Nestor et al. (2011) 
we have carried out the stacking o: 
all the 45 sources in 'Fable (TJ regardless of the contamina¬ 
tion by possible foreground sources. In this case we have a 
clear detection in the U band with S/N = II, oiU = 28.87. 
This translates into an escape fractio n of /J()( = 21%, in 
agreement with previous findings of Iwata et al. (2009)] 
Nestor et al. (20II)| Mostardi et al. (2013) 'The resulting 
ionizing emissivity of L > 0.5L* galaxies (or equivalently 
R < 25.5, our limit) is r_i 2 = 1.24, thus all the LyC 
photons measured in the Lyman forest should be provided 
by relatively bright galaxies, with the consequence of a 
negligible contribution by fainter galaxies or by AGNs. 
Taken at face value, this would imply that only galaxies 
brighter than ~ 0.5L* are contributing to the ionizing 
radiation at z ~ 3.3, thus excluding fainter galaxies and the 
QSO/AGN population as sources of ion izing photons. This 
is at variance with the conclusions of Lusso et al. (2015) 


and Glikman et al. (2011) who found that bright QSOs at 


z ^ 2.4 and at z > 4, respectively, should provide at least 
^ 50% of the HI ionizing photons. _ 


If we integrate the Reddy & Steidel (2009) LF down to 
Mi 5 oo = —17.2, corresponding to a limit of R = 28.5, and 
compute the escaping emissivity at 900 A rest-frame as¬ 
suming = 21%, we end up with r_i 2 = 3.94, well 
above the measured value at z = 3.3 shown in Figl?) We 
can thus conclude that a large fraction of the galaxies with 
detected LyC emission in the VUDS-LBC/COSMOS field 
are plausibly contaminated by for eground interlopers , as 
discussed above and al so shown by Siana et al. (2015) and 
Mostardi et al. (2015) In this case, star-forming galaxies 
alone will probably not be able to keep the Universe ion¬ 
ized at z ^ 3.3. The situation can be different at higher 

redshifts. _ 

Recent evaluations by Finkelstein et al. (2014) show a 
steepening of the faint end of the Lh' with a slope as steep 
as -2 at z ~ 6—7, i.e. close to the epoch of reionization. Very 
high redshift star forming galaxies down to M 1500 = —17 
at z 7 could ionize the I GM if they had an esc ape frac¬ 
tion of ^ 13%. Similarly, Bouwens et al. (2015) modeled 
the redshift evolution of the U V background integrating the 
LF of star-forming galaxies down to M 1500 ~ —13 and as¬ 
suming an escape fraction of 0.10-0.40. These parametriza¬ 
tion allow them to reproduce the observed trend of the 
ionization rate Nion inferred at z > 6 by a number of in¬ 


dependent obseryations_(seetiieir/^ble 1). Using the same 
approach, IDuncan fc Conselice 20151 concluded that star¬ 
forming galaxies at z > 6 down to M 1500 —15 and with 

fesc ~ 0.16 — 0.42 could complete and maintain the reion¬ 
ization process. Given our limits (//// < 2%) at z ^ 3.3, 
an increase in the escape fraction should be invoked for 
star-forming galaxies with lower luminosities and/or in¬ 
creasing redshifts to make these sources the major con¬ 
tributors to the reionization at z > 6 . The observational 
evidence that galaxies at higher redshifts and fainter lu¬ 
minosities are bluer and less dusty than bright and low-z 
galaxies (|Duncan fc Gonselice 20 IS]! and the recent finding 
that ionization para meter can be higher for z > 6 galaxies 
(Stark et al. (2015)), may indicate that high-z galaxies can 
be more efficient ionizing agents than galaxies at z < 3. 

The alternative possibility that the AGN population can 
provide most of the ionizing flux in the Universe cannot be 
excluded, however. Faint multiwavelength surveys based on 
very deep near infrared data and at the same time ultra¬ 
deep X-ray images are discovering a population of faint 
AGNs at very high redshifts (z ~ 4 — 7) that could pro¬ 
vide the required ionizing photons (e.g. Fiore et al. 2012, 
Giallongo et al. 2012, 2015). The main uncertainty even 
for this population relies on the escape fraction of ionizing 
photons from their host galaxies. It is not clear whether 
the AGN escape fraction remains high for progressively 
fainter sources, but, at variance with the galaxy popula¬ 
tion, the brighter AGNs and QSOs having - 100% 
are known to ionize th eir neighborhood even at z ^ 6 
( Worseck et al. ( 2014a)|). In th e local Universe, recent re- 
sults by Stevans et al. (2014) seem to indicate that even 
the faint AGN population has an escape fraction close to 
unity. Accurate measurements of the escape fraction of faint 
AGNs at z > 3 will be important to confirm this hypothe¬ 
sis. 


6.5. Comparison with theoretical predictions 

Theoretical predictions of the ionizing photon leakage from 
star-forming galaxies at high-z are far from being settled, 
despite the large effort put in complex numerical simula¬ 
tions, and different conclusions have been obtained by var¬ 
ious authors on this subject in the past. 


Using high resolution simulations, Gnedin et al. (2008) 
derived an absolute escape fraction of 1 — 3% for galax- 
ies with mass M > IO^^Mq and star formation rate 
^ 1 — ^MqIut at 3 < z < 9. This model is in slight 
agreement with our result on the VUDS-LBG/GOSMOS 
region. They predict also that the escape fraction declines 
steeply at lower masses and SFR, while they found large 
variance among different lines of sight within individual 
galaxies. According to their model, a high-z galaxy can be 
more transparent along the galactic poles than through the 
galactic disk. They explain the low values of escape frac¬ 
tions with a small fraction of young stars located just out¬ 
side the edge of the HI disk. In lower mass galaxies, this 
disk is thicker and more extended relative to the distribu¬ 
tion of young stars compared to massive galaxies, causing 
a progressively low er escape fraction. 

Differently from|Gnedin et al. (2008) , the simulations of 
Yajima et al. (2014)| found that //)// remains roughly con- 
stant at a value of 20% (corresponding to //j/ = 52% for 
E{B — V) = 0.1 and a Calzetti extinction law) from z=0 
to z= 10 , and it does not show any clear dependence on 
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Table 3. The HI photo ionization rate for different LF parametrization and integration limits 


rrel 

Jesc 

LF 

a r 

_i2(Mi5oo < —20.2) 

L > 0.5L* 

r_i2(Mi5oo < —17.2) 

L > 0.03L* 

r_i2(Mi5oo < —14.2) 

L > 0.0016L* 

< 2 % 

Reddy & Steidel (2009)1 

-1.73 

< 0.12 

< 0.38 

< 0.93 

< 2 % 

Sawiki & Thompson (2006) 

-1.43 

< 0.10 

< 0.23 

< 0.31 

< 2 % 

IGucciati et al. (2ol2)| 

-1.50 

< 0.10 

< 0.19 

< 0.26 

< 2 % 

Gassata et al. 

-1.78 

- 

< 0.08“ 

- 

21 % 

Reddy & Steidel (2009) 

-1.73 

1.24 

3.94 

9.73 


The upper limits on and r_i 2 are at Icr. (a) 


Cassata et al. (2011) found a star formation rate density of O.OOSMp, /yr/ Mpc^ 
integrating the LAE LF down to 0.04L*. This SEH density is a factor of 5 below the one by Reddy & Steidel (2009) 


the galaxy properties (i.e. luminosity). This is at variance 
with our results discussed in previous sections (/J|j < 2% 
at L > 0.5L*). 

Using ray tracing simulations with spatial resolution 


of 3.8 pc, Ji-hoon et al. (2013) find = 1.1% for 


spiral galaxy in a dark matter halo of ^ 2 x IO^^Mq. 
They do not predict the dependence of the escape 
fraction from galaxy mass or luminosity. Models based 
on hydrodynamical simulations or radiative trans - 
fer calculations dRazoumov &: Sommer-Larsen (2010 


Paardekooper et al. (2013)] imm & Cen (2014 

Wise et al. (2014)1 generally predict an increase o: 


the escape fraction at fainter luminosities or at lower halo 
masses, with a large value of the absolute escape fraction 
(^ 15 — 80%) at Mi 5 oo ~ [—14, —10] and at high redshifts 
{z > 4). 

Only recently, however, the predictions from cos¬ 
mological hydrodynamic simulations start converging on 
a general consensus that the escape fractions is in- 
creasing towards low mas ses/luminosities. For example 

using high-resolution cosmo- 
ations of galaxy formation, find 


Paardekooper et al. (2015) 
logical hydrodynamic simu. 
that only in 30% of the halos with M 200 = IO^Mq ( corre- 
sponding to a stellar mass of ^ 10^Mq; see their Fig. 5, cen¬ 
tral panels) the absolute escape fraction is higher than 1% 
at z > 9, while > 10% at z=6 and/or at lower masses, 
with a strong anti-correlation between escape fraction and 
halo masses. Again, the population of ultra-faint galaxies 
in their model is essential in providing the LyC photons 
needed to sustain reionization. Our VUDS-LBC/COSMOS 
survey is limited to much brighter objects {L > 0.5L*), thus 
we cannot provide meaningful constrain ts to those mod- 
els. Using 2-D hydrodynamic simulations, Roy et al. (2015) 
find /gj/® ~ 5% for Milky Way (MW) galaxies, with a trend 
of increasing the escape fraction for lower mass galaxies 
if esc ^ 10% for galaxies ten times less massive than the 
MW). 

Using h igh resolution num erical simulations (0.1-4 pc) 
at z > 6, |Ma et al. (2015) find a large time variability 


of /ggg® (0.01-20%), with a mean value < 5% over long 
timescales, without any strong dependence on mass or red- 
shift. They also show that models based on “sub-grid” sim¬ 
ulations, which are not able to resolve the star formation 
processes in dense molecular clouds, usually over-predict 
the HI ionizing escape fraction. 

Phenomenological models (|Haardt fc Madau (2012)] 


IKuhlen fc Faucher-Giguere (2012) 

[Alvarez7 Finlator &: Tfenti (2012)] ' [Fontanot et al. (2014)] ) 
explored different parametrization of the escape fraction 
of high-z galaxies to explain the evolution of the UV 
background and the Thomson optical depth measured 


at that time by WMAP. They all assume a low escape 
fraction at z ~ 3, of the order of /“/(f ^ 3 — 5%, but 
they require a strong evolution with luminosity or with 
redshift, in order to fit the higher redshift data. The escape 
fraction of VUDS galaxies at z ~ 3.3 is still compatible 
(< 6% at 3cr level) with these values. More stringent 
limits and the exploration of the escape fraction at lower 
luminosities and/or higher redshifts can thus provide 
significant constraints to these models. 


7. Summary and Conclusions 


We have analysed the es cape fraction of 45 gala xies selected 
in the COSMOS region (jScoville et al. (2007)| with robust 
spectroscopic redshifts between 3.27 < z < 3.40 from the 
ultra-deep VI MOS spectroscopic cat alog provided by the 
VUDS survey (|Le Fevre et al. (2015)]). Deep U and R band 
images obtained by LBC at the pr ime focus of the LBT 
telescope (observations described in Boutsia et al. (2014)), 
reaching Icr magnitude limits of U ~ 29.7 and R ~ 28.2, 
were used to measure the ionizing escape fraction from the 
galaxy sample down to a few percent level. The main results 
are the following: 


— A sub-sample of ten galaxies with R < 25.5 has a mea¬ 
sured U band detection that corresponds to a LyC emis¬ 
sion of //// > 28%. A detailed analysis of their proper¬ 
ties has been performed computing the ’’local” escape 
fraction in the spatial region of the U band where the 
supposed LyC flux was measured. Unphysical escape 
fractions //Jj] ~ 300 — 1000% have been derived for eight 
out of ten objects. The relative fesc of these galaxies re¬ 
mains above 100% even assuming a transparent IGM 
and non standard stellar populations/lMFs (Poplll or 
top heavy IMF). 

— High resolution HST images in the 1814W band suggest 
that the supposed LyC emission comes usually from dis¬ 
tinct interlopers slightly offset from the positions of the 
star-forming galaxies at z ~ 3.3. It turns out that eight 
out of the ten z ~ 3.3 galaxies show detection in the 
LBC U band contaminated by light from foreground 
galaxies at z < 3. Only two galaxies could be true LyC 
emitters, though the significance of their detection in 
the U band is marginal, at 2cr level only, and for one of 
these two the possibility of contamination from an inter¬ 
loper is still open. Further observations of these peculiar 
sources will be fundamental to define their physical na¬ 
ture. 

— We stacked the U and R band images of the remaining 
clean sub-sample of 37 galaxies at z ~ 3.3 (including the 
two marginal LyC detections) to compute their average 
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ionizing escape fraction. The stacks are undetected in 
the U band down to ^ 31.4 magnitudes (at 1 cr). This 
corresponds to an upper limit of < 2% (1 cr) at 
R < 25.5 or equivalently L > 0.5L*{z = 3). 

Considering the stochasticity of the IGM and the patch¬ 
iness of the ISM along different LoS as two sources of 
uncertainties for the derivation of the escape fraction, 
the relatively large number of galaxies analysed here, 
i.e. 37, allows us to constrain between 1.3% and 

2.8% at 1 cr con fidence level. _ 

Integrating the Reddy & Steidel (2009) 1500 A LF at 
z ~ 3 down to 0.5L* we have derived a significant upper 
limit to the ionizing emissivity and photo-ionization rate 
of relatively bright star-forming galaxies r_i 2 < 0.12 
at 68% confidence level (corresponding to r_i 2 < 0.36 
at 99.7% confidence level). This limit is a factor of 
two lower than the recent estimates derived from the 
IGM observed in the Lyman-a forest spectrum of bright 
QSOs at z=3, as shown in FiglT] We can thus conclude 
that galaxies brighter than 0.5L* give a small contribu¬ 
tion to the UV background needed to keep the Universe 
ionized at z ^ 3. 

Our results differ from those obtained at similar 


redshifts by 

Iwata et al. (2009)| Nestor et al. (2013)| 

Mostardi et a 

. (2013) The reason resides mostly on the 


spectroscopic redshifts are available for all the galax¬ 
ies of our sample and high resolution HST images sup¬ 
ported by spectral synthesis models have been used to 
clean our sample from bona-fide low redshift interlop¬ 
ers. Only few objects are currently genuine LyC emit¬ 
ters, and further analysis is important to understand 
thei r nature. Similar conc lusion s have been obtained 
by Vanzella et al. (2012a) and [Vanzella et al. (2015) 


in the GOOUS-South fi eld, by |biana et al. ( 2015)1 in 


the SSA22 area, and by Mostardl et al. (2015) 


the 
low up 

of all the faint interlopers can assess their nature. 


HS1549-I-1933 field. Future NIR spectroscopic fo 


In summary, considering that the expected AGN contri¬ 
bution at z ^ 3 amounts to F_i 2 = 0.37 (Haardt & Madau 
2012), an additional contribution of F_i 2 ~ 0.06 — 0.33 
is require d to match the value of r_i 2 = 0.791°'^g ob¬ 


served by Becker & Bolton (2013) This additional contri¬ 
bution could be provided by a faint {L < 0.5L*) galaxy 
population with an escape fraction of > 6 — 20% (the 
exact value depends on the assumed slope of the galaxy 
LF), which is at least three times larger than our upper 
limit derived from the brighter population. 

At higher redshifts (z = 4—7), where the contribution of 
bright AGNs to the photo-io nization rate is expected t o be 
progressively less important (Haardt & Madau (2012)), the 
escape fraction of star-forming galaxies is expected to be 
larger than the one measured here at z ~ 3. The galaxy LF 
at z > 6 is quite steep (a < — l-7)i compared to the LFs a t 
3 (e.g. Cucciati et al. (2012)[(Finkelstein et al. (2014) I 


and the contribution of faint galaxies (down to Mjjy ^ “14 
with ~ 10 — 20%) is required to provide enough LyC 
photons to effectively contribute to the reionization of the 
Universe. This hypothesis however would imply a rapid 
change of the escape fraction of ionizing radiation from 
star-forming galaxies significantly fainter than L*, which 
is difficult to explain in the current t heoretical scenar - 
ios for galaxy formation and evolution (Ma et al. (2015) I. 


From the observational point of view, instead, there are 
indications that galaxies are bluer at high redshift and 
faint luminosities, thus plausibly implying younger and 
less absorbed stellar populations, with possible larger 
contribution of ionizing photons than the ones present 
at z ~ 3 (Stark et al. (2015) Bouwens et al. (2015) 


IDuncan fc Gonselice 2015|l . A simiiicant contributions to 
the reionization by different populations is however still 
possible, since the escape fraction of galaxies at z = 
4 — 7 is currently a free para meter. Faint AGNs are 
the most plausible competitors ([Giallongo et al. (2012) 


Worseck et al. (2014b)[ Giallongo et al. (2015)) but more 


exotic alternatives (primordial cosmic rays, FopIII stars 
forming high-mass x-ray binaries or Hawking evaporation 
of primordial mini black holes) could also pla y a role 
(Tueros et al. (2014) [Belotsky fc Kirillov (201^ ). 

In the future, different activities can be pursued to 
shed light on the topic of the escape fraction of high- 
z galaxies and their relation with the reionization of 
the Universe. In particular, we are looking for possible 
faint ionizing galaxy populations combining high spatial 
resolution HST images from the GANDELS survey in 
GOODS-North with deep LBG imaging in the U band 
dGrazian et al. in prep.D in order to measure the ionizing 
escape fraction at z ^ 3 in galaxies (and AGNs) fainter 
than those included in the present VUDS-LBC GOSMOS 
field (i.e. R > 25.5). At relatively bright luminosities, 
the present sample can be significantly enlarged by the 
other two fields of the VUDS survey, i.e. the VVDS 2- 
hour region and the EGDFS, which benefit respectively 
from relat ively deep UV images by GEHT, Subaru, an d 
by VLT (iLe Fevre et al. (2015)| |Vanzella et al. (2010b)|). 
Intrinsically faint galaxies {L << L*} instead could be dA 
tected in LyC thanks t o the large magnification s by strong 
lensing, as prop osed by Vanzella et al. (2012b) and already 
carried out by Amorin et al. (2014)[ The CLASH survey 
(IPostman et al. 20121) and the Frontier Field Initiative by 
HST will provide a large sample of highly magnified galaxies 
at high-z, which will open a new window on the reionization 
process. 
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Appendix A: Other examples of contamination by 
foreground galaxies 

We show here the cutouts in the U and R bands of LBC 
and in the I814W band by HST of peculiar sources of the 
VUDS-LBC/COSMOS field. These sources, in the last two 
subsections of Table [U have been detected in LyC, but 
there are evidences of contamination by foreground sources 
closely aligned with the z ^ 3.3 sources, as mainly dis¬ 
cussed in Section 4. In particular, Fig lA.il shows two ob¬ 
jects with “global” /J®/ > 100% and Fig lA.2l shows the 
example of four galaxies with “global” //// < 100%. In all 
the six galaxies however the “local” escape fraction is in 
the range ~ 200 — 1000%, thus indicating contamination 
by foreground sources. 
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Fig. A.l. The cutouts of sources ID=52336 and 55594 in the U band (left) and in the R band (central panels) by the 
LBC instrument. The right cutouts show the same sky area observed by HST in the I814W filter. The positions of the 
galaxies with spectroscopic redshifts z ^ 3.3 are indicated by the green circles. The size of the cutouts is 12 arcsec for 
source ID=52336 and 25 arcsec for ID=55594, respectively. 
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Fig. A.2. The cutouts of sources ID= 19521, 40998, 48734, and 47357 in the U band (left) and in the R band (central 
panels) by the LBC instrument. The right cutouts show the same sky area observed by HST in the I814W filter. The 
positions of the galaxies with spectroscopic redshifts z ~ 3.3 are indicated by the green circles. The size of the cutouts 
is 15, 10, 12, and 10 arcsec for source ID=19521, 40998, 48734, and 47357, respectively. 


